Proposed low-energy model Hamiltonian for the spin-gapped system CuTe205 
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Using first-principles electronic structure calculations based on the N th order muffin tin orbital 
(NMTO)-downfolding technique, we derived the low-energy spin model for CuTe205- Our study 
reveals that this compound is a 2D coupled spin-dimer system with the strongest Cu-Cu interaction 
mediated by two O-Te-0 bridges. We checked the goodness of our model by computing the magnetic 
susceptibility with the Quantum Monte Carlo technique and by comparing it with available experi- 
mental data. We also present magnetization and specific heat results which may be compared with 
future experimental investigations. Our derived model is in disagreement with a recently proposed 
model for this compound [J. Deisenhofer et a/, Phys. Rev. B,74 (2006) 174421]. The situation 
needs to be settled in terms of further experimental investigations. 

PACS numbers: 75.10. Jm,75.30.Et, 71. 20.-b 



I. INTRODUCTION 

Significant amount of effort both experimental and the- 
oretical has been devoted in the last years to the investi- 
gation of the behavior of low-dimensional quantum spin 
systems^. A crucial piece of information needed in the 
process of understanding these systems is the connection 
between the underlying chemical complexity of the com- 
pound and the spin lattice. Often, this relation is not 
obvious from structural considerations and one needs to 
rely on ab initio based calculations as we will show in the 
present work. 

Recently, in an attempt to analyze the effect of lone- 
pair cations like Se 4+ or Te 4+ on the magnetic dimen- 
sionality of Cu 2+ -based systems, the magnetic properties 
of CuTe205 were investigated 2 . CuTe205 is structurally 
a Cu(II)-dimer system separated by Te ions. Magnetic 
susceptibility measurements show a maximum at T max = 
56.5 K with an exponential drop at lower temperatures 
signaling the opening of a spin gap. The behavior at 
high temperatures follows the Curie law with a Curie- 
Weiss temperature of #=-41 K— , what indicates that the 
dominant interactions in this system are antiferromag- 
netic. Electron spin resonance (ESR) observations sug- 
gest though that the structural dimers of CuTe205 do 
not coincide with the magnetic dimers 2 -. Fitting the sus- 
ceptibility data to different models, such as a pure dimer 
model, the alternating spin-chain model and the mod- 
ified Bleaney-Bowers model, show equally good agree- 
ment of the experimental data 2 . A detailed investiga- 
tion of the magnetic exchange paths using the extended 
Huckel tight binding (EHTB) method suggests that (i) 
the strongest interaction is between Cu ions which are 6 th 
nearest neighbors (Je) and is of antiferromagnetic super- 
superexchange (SSE) type mediated by a O-Te-0 bridge 
and (ii) that the next strongest interaction is of antifer- 
romagnetic superexchange (SE) type within the struc- 
tural dimer (Ji), yielding a ratio Ji/J6=0.59 2 -. Based on 
these findings, the authors of Ref. 2 proposed an alter- 



nating spin chain model as the simplest possible model 
for CuTe 2 5 . 

Given the existence of many possible Cu-Cu interac- 
tion paths in this material whose relative importance 
may not be necessarily captured in EHTB study, we per- 
formed ab initio density functional theory (DFT) calcu- 
lations and applied the NMTO-downfolding technique^. 
This technique has proven to be very successful in deriv- 
ing the underlying spin model of a large number of low- 
dimensional quantum spin systems including cuprates 4 ^, 
vanadates^ and titanatesi Our calculations reveal that 
the strongest Cu-Cu interaction in CuTe205 is the one 
between 4 th nearest neighbor mediated by two O-Te-0 
bridges (J4) followed by the interaction mediated by a 
single O-Te-0 bridge (Je). This is in contrast with the 
fact that J6 was found in Ref.d to be the strongest in- 
teraction. We also obtain that the Cu-Cu interaction 
within the structural dimer unit (Ji) is rather weak as 
opposed to the findings of the EHTB study. The under- 
lying spin model for CuTe205 derived out of our cal- 
culations is therefore different from that suggested in 
Ref.Q. We have also computed the magnetic suscep- 
tibility for the proposed model by performing Quan- 
tum Monte Carlo (QMC) simulations (stochastic series 
expansion&^iii-) . Our results show good agreement with 
the experimental observations. In view of the fact that 
the magnetic susceptibility is often found to be an in- 
sensitive quantity to the details of the magnetic struc- 
ture, we have also calculated temperature and magnetic 
field dependent magnetization as well as the specific heat 
as a function of temperature. These results need to be 
tested in terms of further experimental investigations to 
resolve the underlying microscopic model for CuTe205 
completely. 

The paper is organized as follows: in section II and III 
we present respectively the crystal structure and the ab 
initio DFT electronic structure of CuTe2 05. In section 
IV we discuss the effective model Hamiltonian obtained 
with the NMTO downfolding method. QMC results for 
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FIG. 1: (Color online) Building units of CuTe 2 5 . (a) Cu0 6 - 
distorted octahedron, (b) Cu20io-structural dimer unit, (c) 
TeCU-tetrahedra. 

magnetic susceptibility, magnetization and specific heat 
are described in section V and finally in section VI we 
present our conclusions. 



II. CRYSTAL STRUCTURE 

CuTe205 crystallizes in the monoclinic space group 
P21/oii with lattice parameters a=6. 871 A, 6=9.322 
A, c=7.602 A, and ^=109.08°. It is built out of 
Cu06 distorted octahedra (Fig. []Ji), with six inequiva- 
lent oxygen 01, 02, 03, 04, 05, 05' surrounding each 
Cu(II) ion. Each Cu06 octahedron is elongated along 
the 02-05' axis, with distances dc u -05' =2.303 A and 
dcu-02=2.779 A. The Cu-0 distances in the CUO4 plane 
range from dcu-05 = l-948 A to dcu-03=l-969 A. Two 
neighboring Cu06 octahedra share an edge to form a 
CU2O10 structural dimer (Fig. [TJd). The oxygen octahe- 
dra of two Cu(II) ions within a given structural dimer 
are rotated by 180°with respect to each other. 

The structural dimers form a chain-like structure run- 
ning almost parallel to the crystallographic c axis. These 
chains pile along the crystallographic b axis (Fig. [2]). The 
Tel atoms are situated between two successive Cu(II)- 
structural dimer chains, while the Te2 atoms are located 
in between two CU2O10 structural dimers along a given 
chain. The local oxygen environment of the Te atoms 
form a Te04 tetrahedra (Fig. [It). The layers containing 
these chains in the be plane are stacked approximately 
along the crystallographic a axis with hardly any con- 
nection between the layers. 



III. ELECTRONIC STRUCTURE 

In order to analyze the electronic behavior of CuTe2 05 
we carried out DFT calculations within the local den- 
sity approximation by employing both the Wien2k code 
based on the full-potential linearized augmented plane 
wave (LAPW) method 12 and the Stuttgart TBLMTO- 
47 code based on the linear muffin-tin orbital (LMTO) 
methodic. The calculated bandstructures agree well with 
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FIG. 2: (Color online) Crystal structure of CuTe 2 5 . The 
largest balls represent Tel and Te2. Tel and Te2 are shown in 
grey and black colors respectively. Cu atoms are represented 
by medium sized balls, situated at the centre of the distorted 
octahedra. The smallest balls denote the oxygen atoms. 



each other within the allowed error bars of the various ap- 
proximations involved in these two methods. Fig. [3] and 
Fig. 2] show the non-spin polarized band structures and 
density of states (DOS) respectively of CuTe205. The or- 
bital contributions to the valence and conduction bands 
in the band structure and the DOS were determined by 
defining the local reference frame with the local z-axis 
pointing along Cu-02 bond and the local y-axis pointing 
almost parallel to the Cu-05 bond. 

The predominant feature of the band structure is the 
isolated manifold of four bands crossing the Fermi level 
(Ej), formed by Cu-d x 2_ y 2 orbitals corresponding to the 
four Cu atoms in the unit cell, admixed with O-p states. 
These bands are half filled and separated from the low 
lying O-p and non-cL^^ Cu valence bands by a gap 
of about 0.8 eV and from the Te-p- dominated high lying 
conduction bands by a gap of about 2.2 eV, with the zero 
of energy set at the LDA Fermi level. We note that in the 
low energy scale, the LDA results lead to a metallic state. 
Introduction of missing correlation effects beyond LDA, 
for instance with the LDA+U functional, is expected to 
drive the system insulating, as our LDA+U calculations 
corroborated. 

The valence band shows Cu d xy , d yzi d zx and d 3z 2_ r 2 
character dominated bands in the energy range between 
-2.2 eV and -1.2 eV while the O-p-dominated bands are in 
the energy range between -4 eV and -1.2 eV. The contri- 
bution of 02 character in the conduction bands crossing 
the Fermi level is found to be small compared to other 
oxygen because of the large Cu-02 bond length. The 
Tel-p and Te2-p states show a non- negligible contribution 
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FIG. 3: LDA band structure of CuTe2 0s plotted along vari- 
ous symmetry directions of the monoclinic lattice. The zero 
of the energy has been set up at the LDA Fermi energy. The 
dominant orbital contributions in various energy ranges are 
shown in boxes on the right-hand side. The various Cu-d 
characters are shown in the local reference frame as described 
in the text. 



to the bands crossing the Fermi energy, as pictured in the 
inset of Fig. [H and play an important role in mediating 
the Cu-Cu exchange interaction as will be demonstrated 
in what follows. 
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FIG. 4: (Color online) Partial density of states of Cu-d (in 
black full lines), O-p (in red dashed lines) and Te-p (in cyan 
or gray full lines) orbitals, for CuTe205. The inset shows the 
density of states for O-p and Te-p in the energy range close 
to E/, dominated by Cu-d x 2_ y 2 character. 




IV. LOW ENERGY HAMILTONIAN - 
EFFECTIVE MODEL 

A powerful technique to construct a low-energy, tight 
binding (TB) Hamiltonian starting from a LDA band 
structure, is given by the NMTO downfolding method 3 . 
This method derives a low-energy Hamiltonian by an en- 
ergy selective, downfolding process that integrates out 
the high energy degrees of freedom. The low energy 
Hamiltonian is then defined in the basis of effective or- 
bitals constructed via the integration out process. This 
process takes into account the proper renormalization ef- 
fect from the orbitals that are being downfolded. The 
accuracy of such process can be tuned by the choice of 
the number of energy points (N), used for the selection 
of downfolded bands. If the low-energy bands form an 
isolated set of bands, as is the case under discussion, 
the constructed effective orbitals, the NMTOs, span the 
Hilbert space of Wannier functions or in other words, the 
effective orbitals are the Wannier functions correspond- 



FIG. 5: Bands obtained with downfolded Cu-d x 2_ y 2 basis 
(solid lines) compared to full LDA band structure (dashed 
lines). Eo, Ei and E 2 mark the energy points used in NMTO 
calculation. 



ing to the low-energy bands. The real space representa- 
tion of the downfolded Hamiltonian H=^^j(cjcj + h.c) 
in the Wannier function basis gives the various hopping 
integrals tij between the effective orbitals. 

For the present compound we have derived the low 
energy Hamiltonian defined in the basis of the effective 
Cu-d x 2_ y 2 orbitals, by keeping only the d x 2_ y 2 orbital for 
each Cu atom in the unit cell and integrating out all the 
rest. We show the downfolded bands in Fig. [5] in compar- 
ison to the full LDA band structure. With the choice of 
three energy points Eo, Ei and E2, downfolded bands are 
indistinguishable from the Cu-d x 2_ y 2 dominated bands of 
the full LDA calculation. 

The corresponding Wannier function is plotted in 




FIG. 6: (Color online) Effective C\id x 2_ y 2 orbital with lobes 
of opposite signs colored as black and white. The d x 2_ y 2 
orbital is defined with the choice of local reference frame as 
described in the text. 

Fig. [6l The central part has the 3d x 2_ y 2 symmetry with 
the choice of the local coordinate system as stated above, 
while the tails are shaped according to 0-p x /p y . The 
Cu-d x 2_ y 2 orbital forms strong pda antibonds with the 
0-p x /py tails. 0-p x /py tails bend towards the Te2 atom, 
what indicates the importance of hybridization effects 
from the Te cations and enhances the Cu-Cu interaction 
placed at different structural dimers CU2O10. 

Table-I shows the various dominant effective hopping 
integrals Uj with magnitude greater than 1 meV between 
the Cu(II) ions at sites i and j. The notation for the var- 
ious hoppings is shown in Fig. [71 where the subindex of t n 
corresponds to the n th Cu neighbors. The strongest hop- 
ping integral, £4, is found to be between those two Cu(II) 
ions which are placed at different structural dimers and 
the interaction is mediated by two O-Te-0 bridges, ti, 
which denotes the hopping integral between two Cu(II) 
ions situated within the same structural dimer unit, is 
found to be about 1/3 of the strongest hopping integral 
(U). The second strongest hopping term, fe, mediated 
by one O-Te-0 bridge is about 1/2 of t±. Fig. [7|o shows 
the interaction paths in the ab plane, which are weak in 
general and can be neglected. In particular, we mention 
as examples the hopping integrals ts and £7, which are 
approximately 1/10 and 1/25 of the strongest hopping 
term (£4) respectively. In the following we discuss the 
origin of the various dominant interaction paths. 

A. Strongest hopping term -U 

The strongest hopping term, £4, mediated by two O- 
Te-0 bridges is associated to a Cu-O-Te-O-Cu super- 
superexchange (SSE) path generating the spin-spin cou- 
pling J4. The strength of a SSE interaction through an 
exchange path of type Cu-O-L-O-Cu (e.g., L=Te) de- 
pends sensitively on how the O-L-0 linkage orients the 
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FIG. 7: (Color online) Cu-Cu interaction paths t n . The color 
convention is the same as Fig. [2] 



TABLE I: Cu-Cu hopping parameters corresponding to the 
downfolded C\x-d x 2_ y 2 Hamiltonian in NMTO-Wannier func- 
tion basis. Only hopping integrals of strength larger than 
lmeV are listed. 



hopping 


Cu-Cu distances 
in A 


Hopping parameters 
in meV 


ti 


3.18 


38.7 


h 


5.32 


11.0 


U 


5.58 


112.9 


fe 


5.83 


13.7 


te 


6.20 


59.9 


ti 


6.43 


4.9 



two magnetic orbitals (i.e the d x 2_ y 2 orbitals) centered 
at two Cu sites and also on how the tails of the mag- 
netic orbitals, which have contributions of the orbitals of 
the ligand atom L, are oriented with respect to the cen- 
tral part. In Fig. [8] we show the Wannier function plot 
corresponding to £4, where the effective Cu-d^-^-like 
Wannier orbitals are at the Cu sites between which we 
have found the strongest interaction. The O-p^/p^ tails 
bend towards the Te atoms forming O-Te-0 ligand paths 
which are responsible for the strong Cu-Cu bonding. 



B. Second strongest hopping term -fo 

The hopping integral describes the next strong Cu- 
Cu interaction path, which is mediated via one O-Te-0 
bridge and responsible for the Cu-O-Te-O-Cu SSE in- 
teraction generating the spin-spin coupling Jq. Fig. [9] 
shows the Wannier plots of the Cu-d x 2_ y 2 downfolded 
NMTOs. Here the oxygen tails bend towards the inter- 
connecting Te02 unit to provide an interaction pathway 
between the two Cu sites as in the £4 path. However 
the strength of this interaction is expected to be weaker 



FIG. 8: (Color online) Effective orbital corresponding to the 
downfolded NMTOs, placed at two Cu sites situated at two 
different structural dimer units corresponding to the £4 in- 
teraction. Lobes of orbit als placed at different Cu sites are 
colored differently. Lobe colored black (white) at one Cu site 
represents the same sign as that colored magenta (cyan) at 
other Cu site. 




FIG. 9: (Color online) Effective orbital corresponding to the 
downfolded NMTOs, placed at two Cu sites situated at two 
different structural dimer units corresponding to the t6 hop- 
ping term. Color convention is same as in Fig. [8] 



than £4 since there is only one, instead of two O-Te-0 
interaction paths. 



C. Structural intradimer hopping term-^i 

ti corresponds to the intradimer Cu-Cu interaction 
path which is mediated by 05-05' atoms. In Fig. \W\ 
we show the Wannier function plot, where the effective 
C\i-d x 2_ y 2 like Wannier orbitals are placed at the Cu sites 
of the same structural dimer unit. As we stated above, 
each structural dimer unit is made of two edge sharing 
Cu06 distorted octahedra. In the case of the first octahe- 
dron 05 is situated on the basal plane of the octahedron 



FIG. 10: (Color online) Cu-d x2 _ y2 downfolded NMTOs, 
placed at two Cu sites situated within same structural dimer. 
The 02 sites with long Cu-02 bond lengths have been re- 
moved for better view. Color convention is same as in Fig. [8] 

and 05-p^/p^ form pda antibond with the Cu-d cc 2_ 2y 2 
orbital, whereas 05' is situated at the apical position for 
this octahedron. The reverse is true for the second octa- 
hedron. Therefore Cu-d a .2_ y 2 orbitals of two Cu 2+ sites 
placed at the same structural dimer unit are misaligned, 
what is responsible for the weak Cu-Cu intradimer in- 
teraction. As the Cu-05-Cu and Cu-05'-Cu angles turn 
out to be 96.76 degree, the Cu-Cu interaction within the 
structural dimer, which is weak in general, is in the bor- 
derline where a sign change of exchange interaction from 
antiferromagnetic to ferromagnetic may occur—. 

V. SUSCEPTIBILITY, MAGNETIZATION AND 
SPECIFIC HEAT 

The description of the spin model for CuTe205 as ob- 
tained from the NMTO-downfolding technique turns out 
to be that of a system of coupled dimers in a two di- 
mensional (2D) grid (see Fig. [TTj). In order to check the 
goodness of our proposed model, we have calculated the 
magnetic susceptibility as well as magnetization and spe- 
cific heat properties by considering the following spin- 1/2 
Heisenberg model on a Ni x N2 lattice: 

N i 1 N i 1 

H = Ji ^ ^ (S2i,jS2i+l,j) + J4 ^ (S2i+l,jS2i+2,j) 

j=0 2 = j=0 2 = 

El-! ^L_1 
2 1 2 X 

+ Jq ^ ^ [(S2i,2jS2i,2j+l + S2i,2jS2i+2,2j+l) 

j=0 2 = 

+ (S22+l,2jS22+l, 2j-l + S2i+l,2jS2i+3,2j-l)] (1) 

where Ji, J 4 and J6 are the exchange integrals cor- 
responding to the hopping paths t±, £4 and respec- 
tively. Interestingly, this model reduces to the model 
grid that describes the magnetic behavior of CaCuGe206 
when Ji = 0. In that case, the two dimensional model 



has (using the present notation), two critical points at 
J 6 ^-0.9J 4 and J 6 «0.55JA 




JI .T4 .T6 

FIG. 11: The 2D coupled dimer model shown on a 8 x 4 (32 
site) lattice. Periodic boundary conditions have been applied 
on both directions. The thick, thin and dashed lines represent 
the strongest J4, the next strong J6 and the weak structural 
intra-dimer interaction Ji respectively. The site index k is 
given by k — N\ * / + m, where m runs over number of rows 
in the square lattice (m = 0,1, ... Ni), and 1 runs over the 
number of columns in the square lattice (/ = 0, 1, . . . N2) • 

The analysis of model Eq. [1] has been done by 
the Quantum Monte Carlo method (stochastic series 
expansions^) on a 20 x 20 lattice. While the NMTO 
downfolding technique gives us an estimation for hop- 
ping parameters, it does not provide directly values of 
exchange integrals. The exchange coupling, J, can be 
expressed in general as a sum of antiferromagnetic and 
ferromagnetic contributions J=J AFM +J FM . In the limit 
of large correlation, typically valid for Cu based system, 
the antiferromagnetic contribution, J AFM , is related to 
the hopping integral t by the second order perturbation 
relation J AFM = 4t 2 /U, where U is the effective on-site 
Coulomb repulsion. In absence of a satisfactory approach 
for computing J directly, in the following we considered 
the NMTO downfolding inputs to built up the model and 
starting point for relative estimates of various exchange 
interactions. We define the parameters, 
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FIG. 12: Temperature dependence of magnetic susceptibility 
for CuTe2 05- The circles corresponds to experimental data£ 
and the solid line correspond to calculated susceptibility based 
on a 2D-coupled dimer model. 

which measure the ratio of the inter-dimer J6 and struc- 
tural intra-dimer Ji interactions with respect to the ex- 
change interaction which was suggested from the down- 
folding calculations to be the strongest J4. 

The optimal values of ql\ and C12 as well as the strength 
of the primary interaction J 4 and the effective g factor are 
obtained by fitting the QMC results for the susceptibility: 

X th = ((S z - (S z )) 2 ) (3) 

(/iB an d fc# denote the Bohr magneton and the Boltz- 
mann constant respectively) with the experimental sus- 
ceptibility (in [emu/mol]) at intermediate to high tem- 
peratures viaM x — 0.375(g 2 / J)x th - To simulate the 
low temperature region of the susceptibility data we in- 
clude the respective Curie contribution from impurities 
ag ^cw _ c imp /T. The calculated susceptibility in com- 
parison to experimental data is shown in FigfT2l 

The best fit corresponds to the intra-dimer exchange 
integral J4=92.4 K, very close to the value proposed by 
J. Deisenhofer et af for the strongest dimer coupling. 
The optimal value of the g factor =2. 17 was found to be 
slightly larger than the spin only value of g=2, in agree- 
ment with ESR measurement. The optimal values for the 
coupling ratios in Eq. [2] are found to be a\ = 0.27 and 
c*2 = 0.07, rather close to the estimates, 0.28 and 0.11 
respectively, obtained using the second-order perturba- 
tion relationship between exchange interaction (J) and 
the hopping integral (t) given by the NMTO-downfolding 
study. The weak Ji interaction turned out to be of anti- 
ferromagnetic nature giving rise to a positive sign for a 2 • 
With the stochastic series expansion implementation of 
the quantum Monte Carlo method it is possible to simu- 
late quantum spin models in an external field. In Fig. [13l 
we present the computed magnetization as a function of 
temperature M(T) for various magnetic fields' strengths 
and in Fig. [TH we show the comparison of M(T) for the 
model proposed in this work and the alternating chain 
model of Deisenhofer et al? for H = 12.7 Teslas and 
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FIG. 13: Magnetization plotted as a function of temperature 
for the 2D-coupled dimer model of CuTe20s in an applied 
magnetic field of strengths h/J=0.2, 0.5, 1.0 (bottom to top) 
which correspond to H = 12.7, 31.7, 63.4 Teslas. 



FIG. 15: Specific heat plotted as a function of temperature 
for CuTe2 05 for the 2D-coupled dimer model (solid line) and 
the model of Ref. [2| (dashed line) . 




FIG. 14: (Color online) Magnetization plotted as a function 
of temperature for the 2D-coupled dimer model of CuTe2 0s 
and the alternating chain model of Ref. @ for two values of 
the magnetic field H — 12.7 Teslas and H — 31.7 Teslas. 

H = 31.7 Teslas. Two models show distinctly different 
behavior at moderate to high magnetic fields. 

We also calculated the specific heat C V (T) for both 
models and the results are presented in Fig. [151 While 
the overall qualitative shapes of the C v (T) versus temper- 
ature curves for both models are similar, there are impor- 
tant quantitative distinctions which capture the different 
nature of the models. 

Our above computed thermodynamic quantities pro- 
vide a useful framework to test the validity of our pro- 
posed model in terms of further experimental measure- 
ments. 



VI. CONCLUSION 

The analysis of the electronic structure of CuTe205 by 
first-principles NMTO-downfolding calculations as well 
as the calculation and examination of susceptibility data 
by the QMC method leads to a unique description of 
this system as a 2D coupled dimer model. The strongest 
Cu-Cu interaction is between Cu pairs belonging to dif- 
ferent structural dimer units and connected by two O- 
Te-0 bridges. Two additional in-plane interactions of 
about 1/3 and 1/10 of the strongest interaction have been 
found; the latter one being the structural intra-dimer in- 
teraction. This is in disagreement with recent theoreti- 
cal considerations in Ref. [2|, which suggest the CuTe205 
system as an alternating spin chain system with strong 
intra- and inter-dimer coupling. Based on our proposed 
model, we have also calculated the magnetization and 
specific heat which may be compared with new experi- 
mental measurements. We hope that our work will stim- 
ulate further experimental studies. 
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